Next-to-leading-order predictions for t-channel single-top production at hadron 
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We present the predictions at next-to-leading order (NLO) in the strong coupling for the single-top 
cross section in the t channel at the Tevatron and the LHC. Our calculation starts from the 2^3 
Born amplitude gq — * tbq , keeping the 6-quark mass non-zero. A comparison is performed with a 
traditional NLO calculation of this channel based on the 2-^2 Born process with a bottom quark 
in the initial state. In particular, the effect of using kinematic approximations and resumming 
logarithms of the form \og(Q 2 /ml) in the 2 — > 2 process is assessed. Our results show that the 
2^3 calculation is very well behaved and in substantial agreement with the predictions based on 
the 2^2 process. 
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It is a quite remarkable fact that in hadron collisions, 
top quarks can be produced via electroweak interactions 
at a rate comparable with strong production [J 0, [1]. 
Such unique behavior is mainly due to two factors. First, 
a top quark can be produced together with its SU(2)l 
partner, the bottom quark, with a sizable gain in phase 
space cost with respect to a top and anti-top quark pair. 
Second, among the three possible production channels, 
one entails the exchange of a vector boson in the t- 
channel, leading to an enhanced cross section at high 
energies. 

Given the large predicted cross section, evidence for 
single top production has been actively sought and re- 
cently established at the Tevatron 0, Hj and it will play 
an important role in the physics program at the LHC. 
Single-top production offers, for instance, the only effec- 
tive way of extracting direct information on Vtb [6]. In 
fact, at the Tevatron the prospects for the detection and 
then measurement of the electroweak (EW) production 
cross sections have significantly worsened since the first 
theoretical proposals 7j. The main reason for this was 
an underestimate of the impact of large backgrounds such 
as those coming from W+ jet production (both with and 
without heavy flavors) and from the strong production 
of ti [1]. The situation at the LHC, though bound to 
improve thanks to the larger rates expected, will not be 
qualitatively very different. 

The most accurate analyses for single top are based on 
two essential ingredients. The first is an in situ determi- 
nation of the background rates. Predictions from theory 
are in this case not able to match the needed accuracy. 
The second is the systematic exploitation of theoretical 
predictions for the kinematic properties of signal (and 
backgrounds). This information is encoded via sophis- 
ticated analysis techniques (such as those based on ma- 



trix elements, neural networks and others 0, @])- Such 
methods are crucial in building efficient discriminating 
variables to select the Standard Model signal or possibly 
find indications of new physics effects . 

It is therefore clear that the most accurate predictions 
for the signal, both for rates and kinematic distributions, 
are needed as inputs in these analyses. An intense ac- 
tivity in the last fifteen years has led to increasingly- 
sophisticated predictions at NLO accuracy. Calculations 
have progressed from evaluations of total rates [Hi 11 1 , to 
differential distributions [l2l . [T^L including spin correla- 
tions in production and decay [lj, Qjl [ly, [l7| and finally 
to the implementation of the three production channels 
in a fully exclusive Monte Carlo program 18, lit ]. 
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FIG. 1: Diagrams contributing at LO in the 2 — * 2 (a) and 
2 — > 3 (b) approaches. 

All NLO calculations available so far are based on the 
2^2 scattering process, Fig. [Ua), where a 6-quark 
appears in the initial state 2CJ, |21(. The usefulness of 
such an approach, called the five-flavor (5F) scheme, is 
twofold. Firstly, the calculation greatly simplifies (as we 
shall describe in detail later), leading to straightforward 
calculations and compact results. Second, possibly large 
logarithms of the form log (Q 2 /ml) due to initial state 
collincar configurations with g — > bb splitting are con- 
sistently resummed into the 6-quark parton distribution 
functions leading to an improved stability of the pertur- 
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bative expansion. Effects related to the "spectator &", 
such as the presence of a 6-jet and the b mass, Fig. [TJb), 
only enter at NLO. As a result, most of the current cal- 
culations and corresponding Monte Carlo implementa- 
tions 12, 14, 15, I(| 17, 18, 3| do not accurately model 
such effects. An alternative approach is to consider as 
Born the 2 — > 3 scattering process, Fig. [ljb), keeping a 
finite b mass. In this scheme, called the four-flavor (4F) 
scheme, the b quarks do not enter in the QCD evolution 
of the PDF's and of the strong coupling. The calculation 
of the NLO corrections is much more involved due to the 
inclusion of an additional parton in the final state and 
the presence of a further mass scale. However, features 
associated with the kinematic description of the specta- 
tor 6's can be genuinely investigated at NLO accuracy. In 
this Letter we present NLO results in the 4F scheme and 
compare them to those in the 5F scheme. The two ap- 
proaches, being by definition equivalent, would give the 
same results at all orders in the perturbative expansion. 
At fixed (low) order, however, predictions could in prin- 
ciple differ significantly and the question of the range of 
applicability of each approach is raised. 

A very important simplification in the 2^2 calcula- 
tion at NLO is that QCD corrections completely factorize 
in terms of a light and heavy currents: color conservation 
forbids the interference between diagrams (1-loop/Born 
or real/real) where the light quark and the heavy quark 
lines are connected by a gluon. In addition, real/real in- 
terferences between i-channel (Fig. [Hb)) and s-channel 
(Fig. Ufa)) diagrams vanish. A clear separation between 
the two processes is therefore maintained at NLO. 




FIG. 2: (a) Representative s-channel diagram present in gq — > 
tbq' . Interference of these diagrams with those of Fig. QJb) 
vanishes because of color, (b) Representative t- and s-channel 
diagrams for gq — > tbgq': their interference is suppressed by 
1/Nl 

Quite remarkably, QCD corrections to the 2^3 Born 
still display almost a complete factorization. Most of the 
real emission processes can be uniquely associated (i.e., 
in a gauge invariant way) with either the light quark cur- 
rent or the heavy quark current. Interference terms are 
either exactly zero or color-suppressed by 1/N%. We have 
kept the color suppressed interferences between light and 
heavy currents, but not those between s-channel NNLO 
real corrections and the ^-channel real corrections (de- 
picted in Fig. H^b)). Interferences in qq — > tbq'q sub- 
process between i-channel diagrams and those with an 
on-shell W, qq — > tb(W~ — ► qq') vanish in the W zero- 



width limit and have not been included. However, we 
have checked that all the neglected interferences are very 
small and de facto do not hamper a meaningful separa- 
tion of the channels. 

All the analytic calculations presented in this paper 
have been performed with FORM 22j: tree-level and 



loop matrix elements are computed at the helicity ampli- 
tude level and therefore top spin information is avail- 
able. Tree-level matrix elements and the LO results 
for cross sections and distributions have been accurately 
checked with MadGraph/MadEvent [13|. The loop con- 
tributions have been evaluated in both the dimensional 
reduction and four dimensional helicity schemes, follow- 
ing the procedure outlined in Ref. [121] . Tensor inte- 
grals have been decomposed with the help of a reduction 
routine based on the Passarino-Veltman approach 24 1. 



Scalar integrals have been explicitly computed with stan- 
dard methods and compared numerically with those of 
Ref. [25I ]. Gauge invariance, CP, and kinematic symme- 
tries (nib <-> filt) have been extensively used to check 
the consistency of the calculation. Infrared and collincar 
divergences in the integrated real and virtual contribu- 
tions have been cancelled locally through the use of the 
dipole subtraction technique [26[ in its massive formula- 
tion [53| as implemented in MCFM [II . The contribu- 
tions from the dipole counterterms have been generated 
independently by MadDipole [2^] and checked point-by- 
point in phase space. Finally, as the most important 
check of our calculation, we have derived the results for 
e + e~ — > Z/j* — ► bbg at NLO. This calculation can be 
obtained from ours by simply accounting for the differ- 
ence in the weak couplings, by setting n%b = mt, by ig- 
noring the QCD corrections on the light current and by 
setting the W virtuality positive. We found excellent 
agreement with the results of Ref. [3(J, that have been 
obtained in completely different subtraction and regular- 
ization schemes. 

We now present and discuss the results of the NLO 
calculation in the 4F scheme (Born 2^3) and com- 
pare with those of the 5F scheme (Born 2 — > 2), at 
the Tevatron (pp, -Js = 1.96 TeV) and the LHC (pp, 
y/s = 14 TeV). In our studies we assume m t = 172 GeV, 
m b = 4.7 GeV and use the CTEQ6.6 PDF set 0. For 
the 2^3 calculation we pass to the 4F scheme by 
adding suitable finite terms, as explained in Ref. [32j |. 
As an independent check we have verified that results 
obtained with the explicit four-flavor MRST set [33j 
are fully consistent with those obtained in the corre- 
sponding five-flavor MRST set plus the finite terms. In 
order to perform a fair comparison between the two 
schemes, we strictly follow the approach of Refs. 2(J 21 1 
and, cont rary to most of the available MC implementa- 
tions [SElElEElEiEi in the 5F scheme, we 
compute the real diagrams gq — > tbq' with a non-zero b 
mass. We have checked that this has a negligible effect 
on the total cross sections, as well as in the top and light 
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a 1 ' 1 at NLO (solid) and LO (dashed) in pb 
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FIG. 3: Scale dependence of the 2 — * 2 and 2 — > 3 calculations, 
at LO (dashed) and NLO (solid) order. Factorization and 
renormalization scales in the heavy and light quark lines are 
equal to fi. For the LHC only top production is considered, 
the behaviour of the anti-top being very similar. 



jet distributions at the Tevatron [34j and the LHC. On 
the other hand, the distributions of the spectator 6's are 
significantly affected. 

In Fig. [3] we show the cross sections for top produc- 
tion at the Tevatron and the LHC in the two schemes 
as a function of n/m t , where fi is a common renormal- 
ization and factorization scale. The 4F calculation has a 
stronger dependence on the scale than the 5F one, par- 
ticularly at the Tevatron, which simply reflects the fact 
that the 2^3 Born calculation already contains a fac- 
tor of a s . However, we observe that both calculations are 
much more stable under scale variations at NLO than at 
LO. To establish an optimal central value for the scales, 
we have studied separately the scale dependence associ- 
ated with the light and heavy quark lines. As expected, 
most of the overall scale dependence is inherited from 
the heavy quark line. In the 4F scheme it is minimal 
for scales around mt/2 and m t /A for the light and heavy 
quark lines respectively, which therefore sets our central 
scale choice. In the 5F scheme the scale dependence is 
very mild and we simply choose to* for both lines. 

Table |T] shows the predictions for the total cross sec- 
tions in the two schemes, together with their uncertain- 
ties. The scale uncertainties are evaluated by varying the 
renormalization and factorization scales independently 
between /j,q' h /2 < /j,f,r < ^Ho' 11 with 1/2 < Hf/Hr < 2 
and /J, L /fi H constant. We see that the uncertainty in 
the 4F scheme is larger than (similar to) that in the 5F 
scheme at the Tevatron (LHC). The difference between 
the NLO predictions in the two schemes is rather small, 
with uncertainties typically less than 5% in both cases. 



The exception is the 4F calculation at the Tevatron with 
an uncertainty of around 10%, which is however still of 
the same order as the absolute difference with the 5F 
calculation. The small scale uncertainties together with 
quite modest increases of the cross sections from LO to 
NLO provide a clear indication that the perturbative ex- 
pansions are very well behaved. 

In Fig. [5] we compare NLO predictions for the top 
quark and light jet pseudo rapidity rj and transverse mo- 
mentum pt- To define the light jet we used the kx al- 
gorithm and imposed pr > 15 GeV, AR > 0.7. Results 
are presented as a bin-by-bin ratio of the normalized (4F 
and 5F) distributions. For the LHC only top production 
is shown, with the behaviour of the anti-top very similar. 
Although the predictions differ somewhat, the differences 
are typically at the 10% level and always less than 20%. 
Finally, we study the NLO distributions in r\ and pt for 
the spectator b. We find that the fraction of events at 
the Tevatron (LHC) where the b is central and at high-px 
(\r)\ < 2.5, p T > 20 GeV) is 28% (36%) with a very small 
scale dependence. From Fig. [5] we see that the largest ef- 
fects in the shapes are present at the Tevatron, where the 
spectator b tends to be more forward and softer at high 
Pt than in the 5F calculation (where these observables 
are effectively only at LO). 

We have reported on the computation of the NLO 
corrections to the EW production of top and bottom 
quarks through the i-channel exchange of a W boson, 
keeping the mass of the heavy quarks finite. This allows 
a systematic study of the approximations and improve- 
ments associated with the different schemes for treating 
heavy flavors in QCD. Wc find that the 4F calculation 
is well behaved: it displays a 10% (4%) scale uncer- 
tainty and a modest (very small) increase of the cross 
section from LO to NLO at the Tevatron (LHC). It gives 
rates that are slightly smaller than the 5F predictions 
(by about 6%). The two calculations are consistent at 
the Tevatron, where the uncertainty of the 4F calcula- 
tion is similar to their difference and marginally consis- 
tent at the LHC, where the estimated uncertainties are 
much smaller. Such a difference could be interpreted as 
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TABLE I: Inclusive cross sections (in pb) for t-channel single 
top production at the Tevatron and LHC using (CTEQ6L1) 
CTEQ6.6 PDF's for the (LO) NLO predictions and $ = m t 



m t . 



and fiQ = mt/2 (/Uq = iWt/4) as central values 



for the factorization and renormalization scales for the light 
(heavy) line in the 5F and 4F schemes, respectively. The first 
uncertainty comes from scale variations, the second from PDF 
errors. 
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FIG. 4: Shape comparison for the top quark and light jet 
NLO distributions. The bin-by-bin ratio of the normalized 
(4F and 5F) distributions in 77 and pr is shown. 



FIG. 5: Distributions (upper plots) in r\ and pr and shape 
comparison (lower plots) of the 2 — > 3 and 2^2 calculations 
of the spectator b at NLO. 



resulting from the resummation of the log (Q 2 /m%) in the 
5F calculation, or as an indication of the need for even 
higher order corrections. A NNLO prediction in the 5F 
scheme could help in settling this issue. The 4F calcula- 
tion provides reliable predictions for all the relevant dif- 
ferential distributions, in particular offering for the first 
time genuine NLO predictions for the spectator b rapid- 
ity and pt- A detailed comparison with current Monte 
Carlo based predictions and an extension of this study to 
the production of fourth generation quarks, such as t'b, 
tb' and t'b' are left to forthcoming studies. 
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